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The Analogy 



Instruction Pointer : Machine Code :: 
Ribosome : RNA 
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0x01 F5 
0x01 F6 
0x01 F7 
0x01 F8 
0x01 F9 
0x01 FA 
0x01 FB 
0x01 FC 
0x01 FD 
0x01 FE 



SP points to location next item will 
be placed in block 

SP decreases when you put item on stack 

SP increases when you pull item from stack 

A B 

d ; : 

X 

Y I 
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The Analogies, cont. 

Instruction Pointer : Machine Code :: 
Ribosome : RNA 

The ribosome translates mRNA to polypeptides 
(transcription -> RNA- processing of pre- mRNA 
->mRNA translation) 





R. Garrett et al. The 

Ribosome: Structure, 
Function, Antibiotics, and 
Cellular Interactions 
(2000) 



W- — ® — -^^P 



Figure 1 : The subunits of a ribosome. Side and front view. 

(A) Small subunit. (B) Large subunit. (1) Head. (2) Platform. (3) Base. (4) Ridge. (5) Central protuberance. 
(6) Back. (7) Stalk. (8) Front. 



More Analogies 



I) Canonical shell commands: cp, mv, cc, ar, In, Id, gprof, ... 

I I ) Biological functional elements: DNA polymerase, ATP/GTP 
powered pumps, ribosome, signal transduction pathways, 
measure macroscopic gene expression, . . . 




by Professor L 

2600 Magazine J Winter 2003-2004 

The creation of pi asm ids (or other vec- 
tors) with transgenic pay loads is made possi- 
ble by the existence of DNA splicing 
enzymes. Simple laboratory techniques allow 
the extraction of naturally accuning plasm ids 
from bacteria and splicing the DNA for the 
new gene into them. The hardest part is figur- 
ing out which combination of genes to insert 
into a host /In order to get a. desired effect. 
However, those techniques are bcybnd flic 
scope of this introductory articled. For oar 



cool, I would strongly encourage you to take 



able to achieve. In this article, I will teach 
you how to put the firefly genes that are re- 
sponsible for the firefly's glow into Es- 
cherichia coli (E. coli for short), the 
bacterium that lives in your gut. You're going 
to make intestinal bacteria that glow in the 
dark. 



Folding of Firefly Luciferase 



purposes, -we can just buy plasm ids wilh our 
desired genes from a scientific supply jftfr^ 

a , , ones and cotranslational folding. To provide a physical 

All E, COll plaSIllld Willi thC LuClferaSC geiie we have started an in-vifro study of luciferase folding. We 

in it is Called PUC I S-IUXR, and Can be pUr- ^ unfolding intermediates and investigated their role in 

_ . ... versible off-pathway aggregation . 

chased from many places (see sources 

irase in vitro from the completely unfolded state or from 
LWifolding intermediates is exceedingly slow. Addition of 
reticulocyte cell lysate, but not of bacterial cytosolic extracts strongly accelerates the 
renaturation reaction. A current project aims at the identification, purification, and 
biophysical characterisation of the cellular factors involved. 



section, below ). 




Crystal structure: Conti et al, (1996) Structure 4, 287; PDB entry 1LCI 



References: 




' Herbst R, S chafer U, Se elder R (1997) Equilibrium intermediates in the reversible 
unfolding of firefly (Photinus pyralis) luciferase. J. Biol Chem. 272, 7099-7105. 

; Herbst R, Gast K, Se elder R (1998) Folding of firefly (Photinus pyralis) luciferase: aggregation and reactivation of unfolding intermediates. Biochemist?} 
37,6526-6597. 



hACKER Lab vs. Bio Lab 





• DNA sequence 
synthesis 

• Online can buy for 
$.50/ bp, up to 45 
nucleotide length 
fragment. 

• Buy your own 
peptide/ nucleotide 
synthesizer for 
$500-$25K USD. 



Noble Prize 1984 Bruce Merrifield: 
solid phase peptide synthesis 



lines 



DNA Synthesis - Beckman 
Oligo 1000 



Peptide Synthesis - Applied 
Biosystems 43 1A 
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PCR lets you assemble pieces ad 

infinitum 



Sketch: 




□ 
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Applied BioSystems Real-Time 
PCR machine ($25K-$45K) 
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Engineering 



Engineer a protein 
Engineer an organism 

. . .. Why? 



'There is at present no understanding of this 
hacker mindset, the joy in engineering for its 
own sake, in the biological community." 

-Roger Brent (Cell 2000) 
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Oh, engineered organisms 



• Corn 

• Tomatoes 

• Citrus fruit 

• (...) 

• And our friend, the fruit fly, Drosophila 
Melanogaster 

• Celera, Inc. released information on genomic- 
scale engineering, not available at press time 
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Primary Flows of Information and Substance in 

a Cell 




12 \ Other Cells / 21 C3 - Berlin 



Review 




□tein... hunh? 



Hbutl' I Z Proteim arc built up by d mi rift 
acids that are Jinked by peptide boitd* to form 
■i polypeptide chain. M i Schematic digram t>f 
an amino acid, illustrating thi- rsomunclaluie 
used in this book. A central cirrhon .si-mn <r H i 
li attached to an amino group iNHil, 4* 
dlfbaxyl jjioup [COOH), .1 hydrogen atmti <Hh P 
and a side chain 1K1. ib| In a polypeptide chain 
H:il- cafbosiyl giotip of amino acid j-p has fornix 
n \n.y\\iiv hoi id, n> Lin- .1 nil no group 
Of amirtn acid ip -f 1 - One water molecule is 
eliminated in this process. The repeating units, 
which are- called residues, are divided into 
main-t hai 11 atoms and side chains. I he 
main-chain pun. which is identical in all 
residue, contains a eenlraJ C (1 atom attached 
to an NH iiroup, a t"; r =0 group. and an H ..Hoed 
The side chain R, which Is different for 
UiilcrLTii residues, is hound to the C 4d atom. 
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Why engineer proteins? 



• 1) Engineered macromolecules could have 
experimental use as experimental tools, or for 
development and production of therapeutics 

• 2) During the process of said engineering, 
new techniques are developed which expand 
options available to research community as 
whole 

• 3) By approaching macromolecule as 
engineer, better understanding of how native 
molecules function 



(Doyle, Chem & Bio, 1998) 
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I s this how a "hacker" 
approaches a problem? 

• 1) determine what are elemental tools/ components, 
learn to work with them, develop something new 

• 2) design/architecture of systems 

• 3) note however the physics/ chemistry of proteins, 
the Levinthal paradox, and the amount of effort 
spent on protein folding, i.e. "more time to hack" 

Levinthal Paradox (1968)- 

given a peptide group 3 possible conformations of 

bond angles (p and ip, in allowable regions 
given a protein of 150 amino acids 
= 3 150 possible structures ~= 10 68 

Real folding times are 

time or bond rotation 10 iZ s r» i i r\r\r\ 

0.1 - 1000 sec 

10 68 * 10- 12 s = 10 56 sec=10 48 years 
15 Life on earth 3.8 * 10 9 years 21 C3 - Berlin 



Methods for de novo protein 

synthesis 

Two methods: 

TASP: Template- assembled synthetic proteins 

RAFT: Regioselectively addressable functional ized 
templates 

"Small proteins or protein domains that are structurally stable and 
functionally active are especially attractive as models to study 
protein folding and as starting compounds for drug design, but to 
select them is a difficult task. 

• • • 

Advances in protein design and engineering, synthesis strategies, 
and analytical and conformational analysis techniques allowed for 
the successful realization of a number of folding motifs with tailored 
functional properties." 

16 21 C3 - Berlin 

(Tuchscherer, Biopolymers, 1998) 



Adding functional motifs to stable 

structures 



(Tuchscherer, Biopolymers, 1998) 




FIGL'EE - Locked-ia tertiary fo!d& as exTenrion of tbe TASF concept: by applying die 
principle of a molecular kit. individual secondary structure elements &ucfa as helice& : ^-sheet& : 
turns, and loops aie covalently attached via both cbain end* to appropriately fUacftoaalized 
templates. Tbe re&u=±ig midtibridged molecules, e.g.. locked-in 4-belix <^4) and .^--sheet {0+) 
bundle*. 8fic£-{$2fi} or m £ |re complex arrangement^ e.g. : {^2): ^ molecule: with a built- 
hi pathway for folding. 22 ^ 
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Ligand Binding - protein flexibility 

"In this study, we set out to elucidate the 
cause for the discrepancy in affinity of a 
range of serine proteinase inhibitors for 
trypsin variants designed to be structurally 
equivalent to factor Xa." 
(Rauh, J. Mol. Biol., 
2004) 



Def: Ligand 

Any molecule that binds 
specifically to a receptor site 
of another molecule; proteins 
embedded in the membrane 
exposed to extracellular fluid. 

18 




Figure 1. Stereo view showing the alternative binding modes adopted by inhibitor (4) in thi: two structures 
(a) XtfKYr)bT.A4 and (b) X(SSYI)bT.B4. In a, the g Lycirx: sparer hydrogen bunds trtGLy216; the busy I group of the inhibi- 
tor occupies the S3/S4siie. In b, the glycine spacer hydrogen bunds to Qy2l9; t±ie tiasyL group points away from the 
enzyme, making contacts with a symmetry- related molecule in the crystal (not shown). 



One way to test for ligand 

binding 



(Doyle, Biochemical 
and Biophysical 
Research Comm., 
2003) 

19 



Ligand 



LBD 



LBD 



£3 



Gal4 DBD — » 






I Gal4 DBD 




Gal4 RE 




HIS3 



Fig. 1. Cienetie seleetion using L V. cerevwiae strain PJWMA. In this 
system the nuelear reeep tor's ligand binding domain [LBD) is fused to 
the CiaL4 DKA binding domain (Gal4 DBD). The fusion protein binds 
to the Gal4 response element eon trolling the expression of the HI S3 
gene and if the HI S3 gene is expressed, yeast eel Is are able to grow on 
media minus histidine.. By transforming the expression plasmidseoding 
for the nuelear reeeptors into the yeast strain and plating them onto 
plates minus histidine but eontaining the appropriate ligand, the nu- 
elear reeeptor aetivates transeription of HI S3 gene. In a process 
analogous to elassieal genetie eomplementation, the small moleeule 
eomplements the histidine autotroph, allowing the yeast to survive 
through a process termed lL ehemieal e om pleme nta ti on . / ? 
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Bioinformatics Databases 



Completely sequenced genomes 



Pfam 



SwissProt 



SMART 




COG - Clusters of orthologous groups 



BLAST with CD »-on 
(Conserved Domain) 



PSI-Blast searches the Non-redundant (NR) database 

21 C3 - Berlin 



How to Access the Human Genome (and other 

sequenced genomes) 

ftp: //ftp. ncbi . ni h. gov 



Index of ftp://ftp + ncbi + nih»gov/genbank/genomes 



Up to higher level directory 
^~|A thaliana 
~| Anopheles gambiae 
^Bacteria 
"~|c elegans 
~\D melanogaster 
^~|H sapiens 

IjLeptospira interrogans serovar Copenhageni 
[ | MITOCHONDRIA 
| J>J raus cuius 
^~|F falciparum 
U Plasmodium falciparum 
1 README 



10/17/2003 0:00:00 

5/7/2002 0:00:00 

4/7/2004 18:28:00 

6/14/2002 0:00:00 

10/19/2000 0:00:00 

4/15/2004 0:23:00 

3/22/2004 17:40:00 

11/2/1999 0:00:00 

5/12/2002 0:00:00 

5/11/1999 0:00:00 

10/11/2002 0:00:00 



hs_phsO . f na . gz Survey sequence (appi 
hs_phsl . f na . gz Unordered contigs (ea 
hs_phs2 . f na . gz Ordered contigs (each 
hs_phs3 . f na . gz Finished sequence 



README OLD 



Index of ftp://ftp*ncbi*nih*gov/genbank/genomes/H_sapiens 



"~|R norvegicus 
^~[s cerevisiae 



Up to higher level directory 




1 README 






2 


KB 


10/17/2003 





00 


00 


~| hs phaseO 


f na. 


gz 


94408 


KB 


4/14/2004 


21 


36 


00 


~| hs phase 1 


f na. 


gz 


606370 


KB 


4/14/2004 


22 


35 


00 


~| hs phase2 


f na. 


gz 


48487 


KB 


4/14/2004 


22 


39 


00 


~| hs phase3 


f na. 


gz 


1090520 


KB 


4/15/2004 
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How to analyze a genome, or 
subsequence (pi) 

• 1 st Step: a) Working with unknown protein sequence; BlastP with 
CD on; you're finding similarity to other proteins, similarity of 
entire AA sequence 

b) COGnitor, precomputed BLASTs; metabolic pathways 
annotated; COGnitor more sensitive since 1) found similarities in 
BLAST, pulled them out 2) works on domain level 

• 2 nd Step: SEG (filtering of low- complexity segments); run COILS 
find a- helices; run Signal P find signal peptides; intrinsic properties 
of SMART, DAS 

• 3 rd Step: run PSI -BLAST to convergence; Pfam picks up 60% of 
known homologs (genes with common ancestor); started with few 
genomes 



22 
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How to analyze a genome, or 
subsequence (p2) 

4 th Step: take result from PSI-BL^ST; run Multiple Alignment 
on that; run Consensus 

( http: //www. accel rys. com/ i nsiqht/ consensus, html ) to find 
conserved regions 




JPRED L 

7ilu 



5 th Step: Predict secondary structure: g 
http: //www, compbio. dundee. ac. uk/ —www- i pred/ 

- Prediction method: "J net; two fully connected, 3 layer, neural 
networks, the first with a sliding window of 17 residues predicting 
the propensity of coil, helix or sheet at each position in a 
sequence. The second network receives this output and uses a 
sliding window of 19 residues to further refine the prediction at 
each position." 

- Determine if protein of unknown function, make inferences based 
on structure prediction Q1 rQ 
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PSI -BLAST 



http://www.ncbi.nlm.nih.gov/BLAST/ 

• A normal BLASTP (protein- protein) run is 
performed. 

• A position- dependent matrix is built using the 
most significant matches to the database. 

• The search is rerun using this profile. 

• The cycle may be repeated until convergence. 

• The result is a 'matrix' tailored to the query. 



24 
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Evolutionary Genomics 

• From a phylogenetic tree can infer 
inheritance of proteins, and thereby 
organisms (conserved vs. non-conserved 
domains, etc). 



Definitions: 

homoloqs : if two genes/proteins share a common evolutionary history 
(not nec. same function) 

analogs : proteins that are not homologs, but perform similar function 
paraloqs : products of gene duplication 

ortholoqs : genes that are derived vertically, no guarantee that perform 
25 same function aira-Hnriin 



Three 



| The Three Types of Trees 



♦ Cladogram 

■ Relative recency of 
common ancestry 

■ No measurement of 
time or change 




((A,B),(C,(D,E))) 



The Three Types of Trees 



♦ UltrametriQ tree 



Relative recency of 
common ancestry 

Depicts evolutionary 
time directly as years 
or indirectly as 
amount of sequence 
divergence via 
molecular clock 



time 



A 

■ B 
• C 
•D 



((A,B),(C,(D,E))) 



bs of trees 



| The Three Types of Trees 



♦ Additive tree 



Relative recency of 
common ancestry 

Branch length 
contains additional 
information, typically 
related to the amount 
of change between 
sequences 



i) 



((A,B),(C,(D,E))) 



Rooted Ys IJnroDted Trees 



♦ Clacl(>grains and 
phjdogmm^ can be either 
rooted or unrooted 

■ Cannot define ancestors and 
descendants in the same 
manner 

■ Can still distinguish clusters 

• Particularly useful in looking at 
different functions of related 
proteins. 

■ Root of a tree is not 
necessarily assigned correctly 
by the program. 




Tools that are neat 



BLAST - does the stuff you'd expect it to 

- It finds stuff. 

- There's some math about why that's good, 
it isn't interesting (unless you're a 
statistician, you aren't a statistician, h 
right?). 

- 1 1 works, don't mess with it. 

3DPSSM ^^|H 

- What's a PSSM? ■ mM 



- Whoa, 3D! 

- Does it really work? 



http://www.sbg.bio.ic.ac.uk/~3dpssm/ 

^9 Imperial College of Science. Technology & fTkedidne 

^7X^fBBAif// Fold Recognition Server 



Lawrence Kelle' 



Fold Library Last Updated: Tue Apr 13 06w 



562] Structui 



Disclaimer and Terms of Use 
Last updated: Men, 05 Apr 2004 03:20:53 GMT Visitors To 

Welcome to the 3D-PSSM Web S 



A Fast, Web-based Method for Protein Fold Recognition using 1D and 
3D Sequence Profiles coupled with Secondary Structure and 
Solvation Potential Information. 



Trans- membrane proteins 

- 20AA a- helix and you got a 
transmembrane prot. 

- (see next slide) 



Click on' in the menu bar to the left to 

submit your sequence 
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Identify trans- membrane proteins 







htt p : // www. cbs. dtu . d k/ servi ces/ Si g na I PI 



Figure 12,1 Four different ways in which 
protein molecules may be bound to a 
membrane. Membrane-bound regions are 
green and regions outside the membrane are 
red. Alpha-helices are drawn as cylinders and 
P strands as arrows. From left to right are 
(a) a protein whose polypeptide chain traverses 
the membrane once as an a helix, (b) a protein 
that forms several transmembrane cx helices 
connected by hydrophilic loop regions, 
(c) a protein with several p strands that form a 
channel through the membrane, and (d) a 
protein that is anchored to the membrane 
by one a helix parallel to the plane of the 
membrane. 
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Nobel Prize for Signal Peptides: The 1999 Nobel 
Prize in Physiology or Medicine has been awarded to 
Gunter Blobel for the discovery that "proteins have 
intrinsic signals that govern their transport and 
localization in the cell." The first such signal to be 
discovered was the secretory signal peptide, which is 
the signal predicted by SignalP. 

28 



Three Case Studies 



• Elite Organisms: 

- Single nucleotide change causes measurable phenotypic 
change (i.e. a fish can see different wavelengths of light), 
(Yokoyama eta/. 2000, PNAS) 

• Engineered Biocatalyst Proteins: 

- Diversa Corp, develops methods for high-throughput 
biocatalyst "discovery and optimization" (Robertson eta/. 
2004, Current Opinion in Chemical Biology) 

• Two protein drugs (FDA approved): 

- TPA - Tissue Plasminogen Activator (Genentech 1986) 

- CSF - Colony Stimulating Factor (Amgen 1987) 



29 
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Diversa Corp and High-throughput 


"Biocatalytic technologies will ultimately gain 
universal acceptance when enzymes are 
perceived to be robust, specific and inexpensive 
(i.e. process compatible). Genomics- based gene 
discovery from novel biotopes and the broad use 
of technologies for accelerated laboratory 
evolution promise to revolutionize industrial 
catalysis by providing highly selective, robust 
enzymes." (Robertson et al. 2004, Curr. Op. in 
Chem. Bio.) 
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x Technology 




GigaMatrix™ Automated 
Detection and Hit 
Recovery System 
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Directed Mutagenesis, Enzyme Family Classification by 
Support Vector Machines, and Support Vector Machines 

SVMs 



given here. SVM is based on the structural risk minimiza- 
tion (SRMj principle from statistical learning theory. 2,1 In 
line ally separable cases, SVM constructs a liyperplane 
which separates two different groups of feature vectors 
with a maximum margin. A feature vector is represented 
by x h with physicochemical descriptors of a protein as its 
components. The hyperplane is constructed by finding 
another vector w and a parameter b that minimizes ||w|| 2 
and satisfies the following conditions: 



w ■ x# + b is +1, fory* = + 1 



W ■ X 3 ; + b ± - 1 1 for = - 1 



Group 1 (positive) 

(1) 

Group 2 (negative) 

(2) 



where y 4 - is the group index, w is a vector normal to the 
hyperplane, |fc|/]|w|| is the perpendicular distance from the 
hyperplane to the origin and ||w|| 2 is the Euclidean norm of 
w. After the determination of w and h, a given vector x # 
can be classified by: 



sign[(w - x') + b] 



(Cai, Proteins, 2004) 

Vapnick, V. (1995) 

The Nature of Satistical Learning 
ThetDry. Springer, New York. 



(3) 



In nonlinearly separable cases, SVM maps the input 
variable iiito a high dimensional feature space using a 
kernel function Klx it Xj). An example of a kernel function is 
the Gaussian kernel which has ]>eeii extensively vised in 
different studies: ".su-sw^si 



(4) 



Based on earlier study 2 '" 38 and our own analysis, Gauss- 
ian kernel function seems to produce better results than 
other kernel functions. Linear support vector machine is 
applied to this feature space and then the decision function 
is given by: 



fix) = sign^ 2«2y*Kfo **) + b 



(5) 



where the coefficients and b are determined by maximiz- 
ing the following Langrangian expression: 



i 1 1 



(6) 



i - 1 t - y - 1 

under conditions: 



oti ^ and = 



(7) 



E- 1 



A positive or negative value from Eq. (3) or Eq. (5) 
indicates that- tlie vector x belongs to the positive or 
negative group respectively. 



Legal Problems with BioTech 

Why this is a huge enterprise 



• Approaches to drug 
patenting: 

- Composition of 
Matter 

- Process Patent (i.e. 
especially with FDA 
approval) 

- Structure 
Characterization 

- Use Patent 



• FDA Approval 

- Takes years and 
years 

- A main reason why it 
takes so long for a 
BioTech firms to 
return on investment 
(i.e. target buyouts 
before product) 
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Goals 



• I ntroduce some current issues 

• I ntroduce resources that address some of those issues 

• "I was a teenage genetic engineer" 

- On DNA Polymerase: 

"Because the complexity of polymerization reactions in vitro 
pales in comparison to the enormous complexity of multiple, 
highly integrated DNA transactions in cells, the biggest 
challenge of all may be to use our biochemical 
understanding of replication fidelity to reveal, and perhaps 
even predict, biological effects. I n this regard, any arrogance 
about our current level of understanding should be 
tempered by the realization that the number of template- 
dependent DNA polymerases encoded by the human 
genome may be more than twice that suspected only four 
years ago." (Kunkel and Bebenek, Annu. Rev. Biochem., 
2000) 
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Reading 



THE 

COMMON 
THREAD 



Eugene Koonin: 

• Sequence - Evolution - Function: Computational 
Approaches in Comparative Genomics (2002) 

J ohn Sulston: 

• The Common Thread: A Story of Science, 
Politics, Ethics and the Human Genome (2002) 

Branden & Tooze: 

• Introduction to Protein Structure (1999) 

I ra Winkler: 

• Corporate Espionage (1997) 

• Spies Among Us: The Spies, Hackers, and 
Criminals Who Cost Corporations Billions (2004) 

Presentations from the O'Reilly BioCon 2003: 

$ wget -r -A ppt,pdf 
http://conferences.oreillynet.com/cs/bio2003/view/e_se 
ss/3516 



I OH N HUSTON 
GtORtilNA ftRRV 



Uiiroductfoit to 

Protein Structure 



4* 

t ,fti I j .i„m,J,vi ,'■ |.n'i.>i ffcr.V 




35 



21 C3 - Berlin 



Acknowledgements 



• GIT co-workers: J ohn B, 
Kristin W, Eric D 

• O'Reilly Bioinformatics Con 
2003 

• Some other people. 



Slides: 

http://qiezi.net/ 
email: 

russell@qiezi.net 

36 21 C3 - Berlin 



